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ABSTRACT 


Digital protection of a power transformer is a challenging problem 
because of certain problems associated with transformers. These problems 
include magnetising inrush current, overvoltage inrush current caused by 
switching, tap changing provision, current transformer mismatch etc. 
Traditionally harmonic restraint has been used to overcome the problem of 
inrush currents, whereas some loss of sensitivity has been tolerated to 
provide for the others. 

The choice of filters to be used for implementing harmonic restraint 
is of crucial importance for the reliable operation of the relay. It has 
been shown that good results can be achieved using Finite Impulse 
Response (FIR) filters, and Haar functions have been used to implement 
these. 

By making provision for tap changing and using a variable bias 
factor, the loss of sensitivity has been minimised. 

The feasibility of a drastic reduction in the hardware required for 
on-line, fast and reliable digital protection of a power transformer has 
been demonstrated. The actual hardware implementation uses a single 
chip microcomputer and has only 2 integrated circuits. 



CHAPTER 1 


INTRODUCTION 

Protection of power system components has always been a 
challenging problem due to the inherent requirements of high 
reliability and fast operation of the protective equipment. 
These requirements are often conflicting because high reliabi- 
lity calls for relatively simple hardware, whereas fast opera- 
tion may require more computational speed which necessitates 
complicated hardware. However, with the advent of VLSI 
technology these requirements are being reconciled more and 
more successfully, and exciting new possibilities in the area } 
of protection systems are opening up. 

Highly reliable relays for transmission line protection 
that can operate in a quarter cycle (5 ms) or less have alread r 
been developed using travelling wave phenomenon. The same is : 
not true for generators and transformers, which can be grouped; 
together, because both require unit protection schemes such as 
percentage differential relaying, supplement by other suitable 
protection schemes. Most of the recent efforts to produce 
differential relays have been limited by computational speed 
requirements. In some cases, hardware multiplications and ©the 
means requiring increased hardware have been used to facilitate 
computation, which lead to poorer reliability. 
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The availability of integrated circuits like the 
Intel 2920 Analog Signal Processor should make it possible soon 
to design differential relays on a single chip. This has been 
the chief motivation in inspiring this work, wherein we have 
investigated the possibility of making some compromises to 
realize this aim, and their effect on relay performance. As 
will be pointed out later, digital protection schemes offer 
several advantages over analog ones, and hence the proposed 
scheme is a microprocessor based, on-line digital relay. 

Percentage differential relays for the protection of 
transformers were rep or ted as early as 1931, by R.E. Cordray [l] . 
Harmonic restraint to prevent false tripping due to magnetizing 
inrush current was introduced soon after by Kennedy and Hayward 
in 1938 [2] . A relay incorporating these features and with a 
claimed tripping time of around 20 ms (i.e., nearly one cycle) 
after the occurrence of a fault was reported by Hayward in 
1941 [3] . Following this, several other schemes using analog 
components were added to the literature on this subject [4-9] , 
until G.D. Rockfeller came out with his landmark paper [ 10] in 
1969, discussing the use of digital computer for, amongst 
various other things, transformer protection* This was followed 
by another paper in 1972 [ll] by Sykes and Morrison, proposing a 
concrete scheme which could be implemented digitally and giving 
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simulation results to evaluate its performance. This scheme 
used harmonic restraint to avoid tripping due to the magnetizin 
inrush currents and used infinite impulse response filters to 
separate harmonics. The resulting scheme was slow in operation 
and carried with it all the drawbacks of the analog filters 
mimicked by it. Later attempts recognized this basic limita- 
tion and used digital filtering techniques as they emerged for 
various other purposes like digital signal processing, image 
processing, etc. In 1976, Malik, Dash and Hope proposed the us« 
of Fourier techniques [12] for filtering purposes. Schweitzer, 
Larson and Flechsig used Finite Impulse Response filters [13] 
for this purpose. Several other papers by Degens [14], Thorp 
and Phadke L15], Rahman and Dash [16], discuss new ideas like 
rectangular transforms to minimize multiplications or computa- 
tional effort, but still require complex circuitry for implemen- 
tation. One of the latest papers on the subject, explores the 
use of Haar functions to achieve harmonic restraint [17] and 
proves the feasibility of such a scheme. 

The study of magnetising inrush current and efficient 
algorithms for simulating it are essential for reliable relay 
design. The first detailed investigation of inrush current and 
the factors on which it depends was carried out by Blume et al* 
[19]. Detailed analysis and experiments were carried out to 
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determine the effect of inrush current on relaying [20]. 

Specht suggested some empirical expressions to simulate inrush 
current in a single phase transformer, in 1951 [21] • Sonnemann 
et al. [22] pointed out that the inrush phenomenon in three 
phase transformers was more complicated than single phase trans- 
formers, and some additional considerations were involved. They 
investigated the suitability of harmonic restraint for relaying . 
and determined the worst possible waveforms against which any 
differential relay was to be tested. Several papers on inrush 
phenomenon in three phase transformers followed this work, 
including one by Specht [23]. The most recent approaches reject 
any functional expression for inrush current, and emphasize that 
a set of differential equations has to be solved numerically to 
simulate inrush current accurately [24] . 

References [25] and [26] have been used for infinite im- 
pulse response filter design in Chapter 3. References [27-30] 
pertain to digital signal processing and Haar function filtering. 
The actual hardware design has been carried out using [31] and 
[32]. 

The present work is organized into five chapters, including 
the present one (introduction) • 

In Chapter 2, the basic scheme of percentage differential 
protection has been introduced. The various faults possible in 
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a power transformer have been considered, and the scope of unit 
protection has been defined clearly with respect to these faults. 
Some problems that arise when a percentage differential relaying 
scheme is used to protect power transformers, like magnetizing 
and over-voltage inrush currents, tap changing, etc., have been 
discussed. 

Chapter 3 is devoted entirely to the selection and design 
of the filters that are to be used for obtaining the fundamental 
component of differential and mean through currents, and the 
second harmonic component of the differential current. The 
criterion for selecting any particular filtering scheme has been 
to minimize settling time for a step input. Haar-f unction based 
finite impulse response filters have been chosen finally, and 
their performance has been examined in detail towards the end of 
the chapter. 

The actual scheme implemented is described in detail, step- 
by-step, in Chapter 4. The hardware design and software used are 
also given therein. A good part of the chapter is devoted to 
evaluation of the proposed relaying scheme through simulation and 
actual testing of the hardware. 

Finally, Chapter 5 summarizes whatever we have been able to 
achieve and also what we could not achieve. Our endeavour in 
this work was not to develop a commercial prototype, but only to 
demonstrate the feasibility of a drastic reduction in hardware. 
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while retaining the benefits of earlier designs, and perhaps 
improving upon them. We have succeeded in doing this, and the 
relay developed has been designed on a single chip. However, 
a lot more remains to be done, and some suggestions for future 
research have been given at the end of Chapter 5. The report 
concludes with a list of journals and books we have referred to, 
in the process of developing the proposed relay. 
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CHAPTER 2 

BASIC CONCEPTS OF TRANSFORMER PROTECTION 


2.1 Introductory remarks 

This chapter starts with a brief description of the 
faults that are likely to occur in a power transformer. 
Protection schemes used for protecting the transformer against 
these faults* and also other abnormal conditions such as inrush 
magnetising current (which flows in the transformer when it is 
energised with the secondary open circuited), overexcitation 
inrush current (caused by overvoltages and switching) etc. have 
been outlined. Finally, the chapter concludes with the intro- 
duction of digital protection schemes, on-line as well as off- 
line, which are being explored nowadays. 

2.2 Faults in a power transformer 

Before attempting to discuss the protection schemes used 
for transformers, it is essential to enumerate the possible 
faults which can be encountered in a transformer. According 
to Warrington £18], all possible cases of abnormal operation 
in a transformer can be divided into two categories : 
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(i) External faults : This category includes problems like 
short duration overloading of the transformer, short circuits 
external to the transformer, etc. The primary transformer 
protection schemes should not operate in such cases. Slow 
blowing fuses and thermal relays are examples of the pro- 
tection used against such events. 

(ii) Internal faults : Faults internal to the transformer 
fall in two large groups, namely, 

(a) Faults which can be detected by unbalance of voltage or 
current at the transformer terminals. These faults are 
of a serious nature and can cause immediate damage. 

Some examples of such faults are : 

1) Phase to phase or phase to earth fault on the HV and 
LV (external) terminals. 

2) Phase to phase or phase to earth fault in HV or LV 
windings. 

3) Short circuit between turns of HV or LV winding. 

4) Earth fault on a tertiary winding or short circuit 
between turns of a tertiary winding. 

(b) Incipient faults : These faults are initially of a minor 
nature, and cause slowly developing damage. These faults 
include 
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1) Poor electrical connection of conductors or a core fault 
causing limited arcing inside oil* 

2) Coolant failure. 

3) Inproper oil flow, allowing local hot spots to develop. 

4) Regulator faults and inproper load sharing between parallel 
transformers, causing overheating. 

Some of these faults can be symbolically represented as 
follows : 


HV T£RT. L V 



Fig. 2.1 Transformer faults 

1) Earth fault on HV external connections. 

2) Phase to phase fault on HV external connections. 
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3) Internal earth fault on HV windings. 

4) Internal phase to phase fault on HV windings. 

5) Short circuit between turns of HV windings. 

6) Earth fault on LV external connections. 

7) Phase to phase fault on LV external connections. 

8) Internal earth fault on LV windings. 

9) Internal phase to phase fault on LV windings. 

10) Short circuit between turns of LV windings. 

11) Earth fault on tertiary windings. 

12) Short circuit between turns of tertiary windings. 

13) Sustained system earth fault. 

14) Sustained system phase to phase fault. 

2.3 Protection schemes used 

Group * a' faults require fast disconnection. Relays 
designed for group 'a 1 faults will generally not operate for 
group * b* faults, and relays designed for group ' b’ will 
generally not be fast enough for group 'a* faults. Thus, 
protection schemes against group * a* and group * b* faults are 
not alternatives, but are supplements to each other. 

Past experience has shown that gas actuated relays are 
adequate for group * b* faults. The most popular gas actuated 
relay is the Buccholz relay, which has been used widely for 
several decades now. 
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For group ’a* faults, the use of differential protection 
was suggested as early as, and even before, 1931 [l]. The 
general idea involved in differential protection is to compare 
the operating quantity, which may be voltage or current, at the 
input and output of the equipment being protected* This com- 
parison is a vector conparison, since the quantity being com- 
pared has a magnitude and phase associated with it. Tripping 
is initiated if the magnitude of the difference exceeds a pre- 
determined constant value (or base value). 

For a synchronous generator, this scheme can be implemented 
as shown below. 



Fig. 2.2 Unit protection scheme of a synchronous 
generator 

Both the CTi in the above scheme are identical, having 
the same turn ratio, and 1^ and represent the currents at 
the two ends of the winding. 1^ and Ig are the scaled down 
versions of and Irrespectively. Tripping occurs if 
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where K is the constant base value* 

There are certain problems in this scheme. At large 
currents due to close-by external faults, even though I£ and 
are equal, 1^ and I 2 may differ significantly due to CT 
errors. This has led to the adoption of a modified scheme 
called percentage differential protection or biased differen- 
tial protection, in which the magnitudes of (T^-T 2 ) and the 
mean through current (Ij+T 2 ) are compared. The tripping 
criterion now becomes 

IVT 2 I 2 s I ^5^ I 

where S is called the bias factor. 

This scheme has been used successfully ever since the 1930s. 

For a transformer, the two current transformers have to be 
of different turn ratios, because Ij^ and are not equal for 
normal operation. If Nj/N 2 is the turns ratio of the power 
transformer being protected, then Ij/l 2 “ R^A^l * if we ignore 
the magnetising current, which is small as compared to load 
currents. As a result, the turn ratios of the CTs are in the 
same proportion as the turns ratio of the power transformer 
being protected. 
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Further, if the winding connections of the power trans- 
former are different on the two sides (e.g., star-delta or 
delta star, which is normally the case) , then there will be a 
phase shift between the primary and secondary current. This 
is taken into account by connecting the CTs in the opposite 
way on the two sides, so that the phase shift is annulled. 

All these points are illustrated in the scheme shown 
below: 



• -sY A 


TRIP IF 

X *1 HhX 0 

I i 1 -r 2 | > s | i 


Fig. 2.3 CT connections for protection of a Y- A 
connected transformer 


Let us examine the threshold characteristic of a percen- 
tage differential relay now. Let us assume that the current 
I 1 is the reference, i.e., = 1^/0°. Further, suppose 

T g = I 2 Zg. Then, 



I x +I 2 



and I x “l2 


= Ii/o + I 2 Z0 

+I 2 cos0 + j I 2 sin0 
* I x -I 2 cos0-jl 2 sin0 


% 


0 



1A 


The threshold condition is 


|Xi-x 2 l =jS ix x +x 2 1 


Squaring both sides, we obtain 

IT T |2 1 J1 »=• . =■ |2 


Hi-12 1 - f sr ir 1 +r 2 r 


Using the expressions obtained earlier, this becomes 


(1^2 cos0)^ + (I 2 sin0) 2 -f-CUj+Ig cos0) z + (l 2 sin0) z ] 


Dividing by 1^ and simplifying, we obtain 

Ip p e2 Io 

U- f-)(j^) 2 + 2(1+ cos0 = f- - 1 

S 2 

Dividing by (1 - 4") » we obtain 


(jS.) 2 + 2(4±4) &) COS0 = -1 

X 1 4-S X 1 


4 , p 

On adding ( • k ) to both sides, this becomes 

4-S 


'2*2 . ^/ X 2n /4 +S 2 n , /4+S 2 \2 / 4+S 2 \2 


+ 2( f i)(^%) cos0 + - 1 * 

X 1 X 1 4-S 2 4-S 2 4-S 2 (4-S 2 ) 2 

4+S 2 4S X 0 

Now defining C = — %, r = — =*- » and p - *=• and substi* 


4— S 


4— S 


tuting these, we get. 
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p 2 -2pC cos0 + C 2 = r 2 


This is the equation of a circle in polar coordinates. 

The bias factor S is generally around 0.05 (5 yi) for 
synchronous generators, whereas it is kept between 0.1 and 0.4 
for transformers. The reasons for using a higher value of S 
in the case of transformers will be discussed briefly now. 

(i) Magnetising current : The primary current contains 
magnetising component which may be up to 5 yi of full load 
current. This current is not there on the secondary side, hence 
a permanent mismatch is created. 

(ii) Current transformer characteristics : Usually it is not 
possible to find current transformers having the f exact turns 
ratio required, since only some standardized models are availa- 
ble. Hence turns ratio requirements can be satisfied only 
approximately, and this leads to some difference in the currents 
being supplied to the relay even under normal conditions. 

Moreover, since the current transformers employed are of 
different turn ratios, their cores may not get saturated by 
the same amount in the event of high mean through current 
resulting from a close by external fault. This may lead to 
some error in the currents being monitored. 
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(iii) Tap changing : Normally power transformers are 
provided with tap changing mechanism while CTs do not have 
such provision. As the transformer turns ratio is changed 
from the nominal value, an unbalance current flows through 
the differential relay. 

(iv) Magnetising inrush current : To understand the phenomenon 
of inrush magnetising current, let us concentrate on a single- 
phase transformer first. 

Under steady state conditions, the flux 0 in the core of 
a transformer with its secondary open circuited, is in qua- 
drature with the supply voltage waveform is as shown below. 



Fig. 2.4 Steady state waveforms of voltage and flux 
in a transformer 

A sinusoidal variation of supply voltage at any instant 
requires a sinusoidal variation of flux in accordance with 
Lenz* s law, i.e.. 
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v = -N 



where N = number of turns, t = time. 

The amplitude of flux variation required for a given amplitude 
of supply voltage can be obtained from the steady state rela- 
tionship 


v 

rms 


4.44 f0 maV N 
max 


As a result, if the core has a residual flux 0„ before 
switching on, and the transformer is now switched on at any 
instant, the resultant flux in the core will be a sinusoidal 
waveshape added to a constant flux 0 r . To clarify this, 
suppose the transformer is switched on at a voltage zero. The 
maximum flux level attained now would by (0 r ±20 maX ) as shown 
below. 
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Since most power transformers operate fairly close to 
saturation flux level, the asymmetrical flux required ini- 
tially on energisation may force the core deep into satura- 
tion as depicted above. This would lead to an abnormally 
high requirement of magnetising current, the exact shape of 
which could be determined from the excitation characteristic. 

In fact, inrush current magnitudes of upto 15 times the rated 
current have been encountered. 

Another problem is created by the fact that this current 
decays very slowly due to the large magnetising inductance of 
the windings, taking several seconds to decay to acceptable 
levels in some cases. 

The phenomenon of inrush current in three phase trans- 
formers is more complicated, because the current flowing in a 
winding is dependent on the currents in other windings also. 

The magnitude of inrush current depends on several factors 
such as the steel used for core construction, the level of 
residual magnetism as well as its polarity, the instant of 
switching in, and so on. Several detailed investigations 
have been carried out on the subject [19-24], and the phenomenon 
has been adequately explained for single phase as well as three 
phase transformers. 
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For relaying purposes, it is sufficient to note that 
the second harmonic content of an inrush current waveform 
is never expected to be less than 16>£ [22] of the fundamental, 
for single as well as three phase transformers. This knowledge 
has been used to discriminate between fault currents and inrush 
currents successfully for a long time and gives satisfactory 
results. 

(v) Over voltage inrush currents : Overvoltages as high as 
20 to 50 •/. result from sudden loss of load. The magnitude of 
such overvoltages is determined by several factors such as 
generator excitation system response, the line length remaining 
connected to a generating station after loss of load, existence 
of faults preceding load rejection, system shunt reactances, 
total generating capacity, etc. Since modern transformers 
operate near saturation flux levels, such overvoltages cause 
a very large increase in excitation current, as large as 10 to 
100 times. This increase may be sustained for long intervals, 
and may cause false operation of a differential relay. 

An increase in magnetizing current means a large increase 
in the third harmonic. This could be used to recognize an over- 
voltage inrush current. Unfortunately, it is not always possible 
to monitor the third harmonic because it gets trapped inside a 
delta connection. The only relaying option left is to monitor 
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the fifth harmonic, which is also present in excitation 
current to a significant extent. Third harmonic monitoring 
would necessitate additional hardware. 

The imbalance between primary and secondary current caused 
by the above factors even in the absence of an internal fault 
necessitates the use of a higher value of bias factor,s,for 
transformers. As hinted above, false operation due to magneti- 
sing inrush current and over-excitation inrush current caused 
by overvoltages/switching is prevented by monitoring the second 
and fifth harmonics of the differential current, respectively. 

As in other fields of protection, initial differential 
relays designed were electro-mechanical. These relays suffered 
from several drawbacks, such as high burden on instrument trans- 
formers, slow operation, contact racing and contact pitting, 
high maintenance requirements, and false operation due to stocks 
or mechanical vibrations caused by external factors. With the 
development of electronic valves, some electro-mechanical 
relays were replaced by electronic ones. Electronic relays 
imposed low burden on instrument transformers, had no contacts 
or moving parts, were easy to maintain and very fast in opera- 
tion, and had none of the drawbacks of electro-mechanical 
relays. Inspite of this, electronic relays did not find very 
easy acceptance, because they too had several drawbacks like 
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high quiescent power con sump tion , very high battery voltage 
requirement, voluminous circuitry, and uncertain life of valves. 

As a result, valves were abandoned as soon as the transis- 
tor and other static components were developed .Subsequent 
experience showed static relays to be stable, reliable, compact 
in size, cheaper, and faster in operation. They can be made 
shock-proof and require very little maintenance and repair. 

They are more sensitive and permit the use of smaller CTs. More 
sophisticated characteristics can be realized. The only major 
drawbacks are complexity and change in transistor characteristic 
with time (aging). 

Both these problems have been overcome with the introdu- 
ction of digital relays. There has been an increasing tendency 
to use these ever since the advent of microprocessors. Initial 
schemes proposed the use of an off-line central digital computer 
for protection. More recent attempts focus on inexpensive, 
reliable, very fast on-line microprocessor based relays. The 
fact that different relay characteristics can be realized with 
the same hardware (or with a minimal change in hardware) adds 
to the charm of these relays. Moreover, fault conditions could 
be simulated and relay performance monitored, making these 
relays self checking, thus reducing maintenance and increasing 
reliability. 
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2.4 Scope of the present work 

Our endeavour has been to design and fabricate an on-line 
percentage differential relay for transformer protection using 
digital signal processing techniques. The Analog Signal 
Processor IC Intel 2920 has been used. To define the scope of 
the work clearly, we must make the following clarifications at 
the outset : 

i) Tertiary windings have been ignored. 

ii) Over-voltage inrush current phenomenon has not been taken 
into account due to limited processing ability. 

iii) Separate relays for detecting external faults, incipient 
faults, and earth faults have been assumed. 

iv) The tripping signal given -is an analog output going high 
when tripping is required. 
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CHAPTER 3 

FILTER DESIGN 


3*1 Introduction 

It was mentioned in the preceding chapter that the 
tripping criterion for a differential relay is given by 

— I,+Y,o 

Uji-M 2 s l-V 2 ! 

Implicit in this statement was the assumption that 1^ 
and I 2 are sinusoidal waveforms, and hence phasors could be 
used to express them. However, in practice, these currents 
contain the fundamental component as well as several harmonics, 
and the resulting waveforms are periodic but non-sinusoidal. 

To implement the above mentioned tripping criterion, we must 
filter out the fundamental components of (Y^-Ig) and (Yj+Yg) • 
Denoting these as i anc * (Yj+I^)^, we can write the 

the tripping criterion as 

i. s 2 | 1 

This expression eaphasizes the fact that we compare the 
rms value of the fundamental components of these currents, which 
is obtained after appropriate filtering. 
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In order to detect the inrush magnetising current, we 
need the second harmonic component (T^-^) as well. On the 
basis of results obtained by Rockfeller et al. [22], we can 
specify an inrush current as being the condition when 

ITj-ijij > 0.16 ITj-Tgli 

where is seconcl harmonic component and (I^-^) , 

is the fundamental component. 

Combining this knowledge with the tripping criterion 
discussed earlier, we can write the final tripping condition as, 

^1 -I 2^1 — S l“TT“ll and ^1~^2^2 - 0,16 

Clearly, such a tripping criterion would give an inrush- 
proof relay, i.e., a relay which would not trip due to magneti- 
sing (transient) inrush current. 

The tasks to be performed can now be visualised as below : 

i) Obtain samples of I ^ and I 2 , convert them to digital form, 

I,+I 2 

and construct (I^-Ig) and ( — • 

X l +I 2 

ii) Extract rms values of (Ij-Ig)!* and ^ — "T* * 

iii) Test if the tripping criterion is satisfied. If yes, 
initiate tripping. 
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The relaying scheme, then, is represented by. the following 
block diagram : 



Fig. 3.1 Block diagram of a percentage differential relay 

In this chapter, we shall concentrate on the design of 
the three filters required in the above block diagram (Fig. 3.1) 
After discussing their specifications, we demonstrate that 
infinite impulse response digital filters cannot give satis- 
factory results. Then we move on to the finite impulse response 
(FIR) digital filters and select the most suitable from amongst 
them, and further explore the chosen filters in some detail. 

3*2 Filter specifications 

It is evident that the specifications for the filters 
extracting (Xj-Ig)^ and (Ij+^^l wil1 ^ visually identical. 
These filters will be referred to as the '50 Hz filters* 
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* 

hereafter. The filter for extracting (I 1 -I 2 ) 2 will be 
referred to as the *100 Hz filter*. 

Let us take up the 100 Hz filter first. This filter 
must reject 50 Hz and 100 Hz frequencies adequately, as these 
are contained in the normal current as well as the fault 
current. How much attenuation is adequate for our purpose? 

To answer this, suppose (I^-I 2 ) a P ure s 3- ne wave of 
frequency 50 Hz, having rms value equal to unity. If the 
100 Hz filter attenuates 50 Hz by 20 dB, then its output will 
be -20 dB, or 0.1, while the 50 Hz filter gives an output of 
1.0. The second harmonic is perceived as being 10/ of the 
fundamental component, even though it is totally absent. To 
ensure reliable tripping, we would definitely prefer better 
Identification of the second harmonic component. 

If the 100 Hz filter attenuates 50 Hz by 30 dB, then 
the * leakthrough* as calculated above will be only around 
3/ f which appears to be reasonably small. Let us, then, 
assume that we require around 30 dB attenuation of 50 Hz by 
the 100 Hz filter. 

Since a large amount of third harmonic (150 Hz) may also 
be present in several cases, let us assume that the same atte- 
nuation is required for 150 Hz also. The 100 Hz filter can, 
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then, be specified as one which has a centre frequency of 
100 Hz, and provides nearly 30 dB attenuation at 50 Hz and 
150 Hz. The exact shape of the frequency response is not of 
any great importance to us, but all harmonics and d.c. must 
be attenuated sufficiently. 

The constraints on the 50 Hz filter will be much less 
stringent, because the 50 Hz component dominates all other 
components of the signal in a power system. Hence, we can 
expect the 100 Hz filter to be the bottleneck in inproving 
response time. Let us explore the possibility of using the 
infinite inpulse response technique to design the 100 Hz 
filter, without worrying about the 50 Hz filter right now. 

3.3 Choice of filtering technique 

Two classes of digital filters are defined, viz: non- 
recursive filters, where the output is a function of only the 
previous and present inputs, and recursive filters, where by 
using feedback, the output becomes a function of past and 
present inputs and outputs. The non-recursive filter generates 
a finite impulse response and is therefore called a Finite 
Impulse Response (FIR) filter. Recursive filters, because of 
feedback, have infinite impulse responses and are referred to 
as Infinite Impulse Response (HR) filters. 
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The traditional approach to the design of HR digital 
filters involves transforming an analog filter into a 
digital filter that meets the prescribed specifications. 

This approach is useful, because the art of analog filter 
design is highly advanced. The conversion of an analog 
filter into a digital filter is acconplished in a number of 
ways; inpulse invariance, direct transformation using s- to 
z-plane transforms, conversion of differential equations to 
difference equations, and direct synthesis. Some popular 
direct transforms are the matched z, left integration, Tra- 
pezoidal, Bilinear, and Prewarped Bilinear Transforms. 

On the other hand, FIR filters offer several advantages. 
They can be designed to have exactly linear phase, thus 
preventing any phase distortion of the input signal. This is , 
important for applications such as data transmission and 
speech processing. Due to their non-recursive structure, FIR 
filters are always stable. However, a strong disadvantage 
is that a large number of delay elements are needed to obtain 
a sharp cutoff, thereby requiring a large amount of .processing. 

Since we had limited processing ability at our disposal, 
we studied the possibility ©f using IIR filters first. Since 
analog filters are geometrically symmetric we can have 30 dB 
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attenuation (for the 100 Hz filter) at 50 Hz and Y"100 2 /50 = 

200 Hz or at 150 Hz and y 100 2 / 150 = 66.67 Hz. The frequencies 
chosen are 66.67 Hz and 150 Hz for the purpose of this study. 
The types of filters investigated by us are Butte rworth, 
Chebyshev, Bessel, and Legendre filters. Some others have 
been left out because they either appeared to have a very low 
cutoff rate, or had very large settling timec (e.g., elliptic 
filters have very large settling timer). The quantity of 
interest to us is the settling time for a step input, because 
we cannot take a tripping decision reliably until the filter 
outputs have settled down. 

Table 3.1 summarizes the results of our investigation of 
FIE filters. From the table, it is clear that the 100 Hz fil- 
ter, if designed using analog techniques cannot have a settling 
time less than around 3 cycles (56.2 ms). 

Clearly, a Legendre filter of the fourth order gives 
minimum settling time. This minimum settling time (56.2 ms) is 
of the order of 3 cycles. If we desire tripping within a cycle 
after the occurrence of a fault, HR filters will, obviously, 
not perform satisfactorily.- Let us focus our attention on FIR 
filters now* 

Traditional FIR filters are digital equivalents of a 
tapped delay line analog filter, or transversal filter. Besides 
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Table 3,1 


Filter type 
Butterworth 


Chebysehev, O.Ol dB 
ripple 


Chebyshev, 0.1 dB 
ripple 


Bessel 


Legendre 


Order 


-3 dB BW 5?< selling time 


2 

3 

4 

5 

6 
7 


15.0 

25.0 

32.0 
37.5 

42.2 

45.2 


97.2 

74.2 
93 .2 
89.8 
88.9 

107.6 


2 

3 

4 

5 

6 

7 

8 


15.0 

25.0 
37.5 

44.1 

50.0 

55.0 

59.0 


56.7 

72.2 

74.1 

72.5 

71.5 
85.0 
85.9 


2 

3 

4 

5 

6 
7 


14.4 

28.2 

38<>9 

47.8 

53.6 

58.1 


103.2 

65.9 

75.1 

69.8 

66.4 

83.5 


2 12.5 

3 18.7 

4 21.4 

5 25.0 

6 


112.3 
82.4 
78.8 
72 .7 
77.3 


2 

3 

4 

5 

6 
7 


13.3 

27.3 

37.7 

45.7 

51.9 

56.9 


62.1 
62.0 
56 .2 
72.0 
69.9 
80.4 


Table 3.1 A conparison of various analog filters with centre 
frequency 100 Hz and 30 dB attenuation points at 
66.67 Hz and 150 Hz. Data obtained from references 
125] and [26] 
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requiring a large number of delay elements (memory locations) 
they require a lot of computational effgrt because of the 
weighing factors involved. Since multiplication is a time 
consuming operation, a simple hardware realization of such 
filters becomes difficult. 

Another option left open is to use any one of the trans- 
form domains employed in image processing. A brief comparison 
of these transforms is given in Table 3.2, which has been 
reproduced partly from [27] . 

A Fourier transform would give the harmonic amplitudes 
directly. If we use any other transform and wish to obtain 
the harmonic contents, some subsequent processing will be 
required. Even so, the table 3.2 indicates that Haar trans- 
form achieves such a drastic reduction in computational effort 
that there may be a reduction in the overall effort required. 

i 

This is borne out by previous attempts using FFT [12] and Haar 
transform L17]. 

In view of the facts mentioned above, it will be desirable 
to use Haar transform for filtering what follows is a brief 
review of Haar functions and their use in digital filtering. 

3 .4 Haar transform-based filters 

The Haar functions form a comple te set of orthogonal 
rectangular functions. They were established by the Hungarian 
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Table 3.2 


Transform Domains 

^ansform name 

Singular value decomposition 

Karhunen-Loeve 

Hotelling, principal 
c oraponents, factor 
analysis 

Cosine 

Fourier 

Slant 

DLB 

Walsh 

Haar 

Hybrid 

DP CM 


for Coding 

Order of computations 
4N 3 

2N 3 

2N 2 log 2 N 

2N 2 log 2 N (complex) 
2N 2 log 2 N 

2N 2 log 2 N (integer) 
2N 2 log 2 N (additions) 
2(N-1) (additions) 
N 2 logN + nN 2 * 

2nN 2 * 


te : For the entries marked with an asterisk (*) , n is proportional 
to the order of predictor used, N is the number of samples foi 
slant, Walsh, Haar? N = 2 k where k is an integer. For cosine, 
Fourier, DLBj order of computations is more than the estimate 
given if N ^ 2* with k being an integer. 
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mathematician Alfred Haar in a paper published in 1910 [28] . 
He described a set of orthogonal functions, each taking 
essentially two values, and yet providing an expansion of a 
given continuous function, which could be made to converge 
rapidly and uniformly* This was a property not obtainable 
by any other set of non-sinusoidal orthogonal functions at 
that time. However, very little practical use was made of 
these functions for over half a century until the 1960s, when 
they were seen to provide some computational advantages in 
certain areas of communication, image coding, and digital 
filtering. 


If we consider a time base 0 < t < t, then we can define 
Haar functions as. 


HAR(2 p +n,t) 






0 


for | p < t < (n+ ^)/2 p 
for (n+ j)/2 p <, t <_ (n+l]/2 p 
elsewhere • 


The first 8 Haar functions are shown in Fig. 3.2. 
It can be shown from the definition that 


1 

/ HAR(m,t) HAR(n,t)dt 
0 


1 for n = m 


1 0 for n m 
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Fig. 3.2 The first 8 Haar functions 

This constitutes the property of orthogonality. The 
proof of completeness is given by Haar [28]. Thus, Haar 
functions as defined above form a complete set of orthogonal 
functions. 

Any given continuous function, periodic with a time 
period T - 1, can be synthesized from a Haar series by 

00 

f(t)= I C HAR(n,t) (3.1) 

n=0 n 

1 

C = / f(t) x HAR(n,t) dt 

n 0 


where 


(3.2) 
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The convergence features of Haar series have been 
examined by Alexits [29]. It has been shown by Shore [30], 
that, a comparatively small number of terms is necessary to 
approximate a waveform using the Haar series. 

If a function f(t) is specified by N samples 
(X^,X 2 » ... , Xj^j) during a time period, we cannot use a Haar 
series for it. However, we can define the discrete Haar 
transform and its inverse from equations (3.1) and (3.2) as. 



, N-l 





= M 1 C i 
N i=0 1 

HAR 

(n. 

i/N) 


N-l 




c i 

= £ X 

n=0 n 

HAR 

(n. 

i/N) 

i,n 

= 0,1,2, 

. • • , 

N- 

-1 


Written in the matrix form, these equations become , 


" ft ^NXN [C] Nxl (3 * 3 > 

and 

tC W ■ [H] n"i tX] N*l (3 ‘ 4) 

where [H] is the Haar function matrix and [H]" 1 is its 
inverse. [H]”* is simply given by the transpose of [H], 
because of the orthogonality property. For N = 8, the matrix 
[H] would be 
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1 1 
1 1 
f2 f2 


[ h] = 


O 0 
2 -2 


O 0 
0 0 
O O 


1111 
1 1-1-1 
f 2 f 2 O O 

0 O V2 f2 

0 0 0 0 

2-200 
0 0 2 -2 

0 0 0 0 


1 1 
-1 -1 
0 0 
V2 1[2 
0 0 
0 0 
0 0 
2 -2 


(3.5) 


2 

If this transformation is carried out directly, N 
operations are required. However, if a factorization 
algorithm similar to the one used for fast Fourier Transform 
is used, the total number of operations can be reduced to 
2(N-1) . A flow diagram for the 8 point fast Haar Transform 
is shown in Fig. 3.3. The solid lines show addition and 
the dotted lines show subtraction at the nodal points. The 
multiplications have all been delayed until the transforma- 
tion is complete. 



= Cc 


= c. 


* c»/Jz - c 3 


« c 4 /z _ c; 


„ / 
-5 




^ f 
*7 


Fig. 3.3 Flow diagram for an 8 point Fast Haar Transform 
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Once the Haar coefficients have been determined using 
the above scheme, the Fourier series coefficients can be 
obtained as derived below : 


Let f(t) = sin ^ t + f2 F 2 cos ~ t + ... 


+ ]f2F^ sin -j t + V"2F^ cos —j t + ... 


We can assume, without any loss of generality, that the 
fundamental time period T * 1. The coefficients F^Fg, ••• 
are given by 


1 

F, - f2 f f(t) sin 2%t dt 
1 0 

1 

F« = f 2 / f(t) cos 2wt dt (3.6) 

z 0 


• • 

• # 

# ♦ 


Now f(t) 
that 


E C HAR (n,t). It follows from equation (3.6) 
n=0 n 


X O® 

F, = V"2 / [ £ C. HAR(n,t)] sin 2wt dt 

A 0 n=0 n 

1 1/2 

* f2 [C / sin 2*t dt + C, .2. / sin 2*t dt + ... ] 
°0 0 


0.9C X + 0.1865 i-C^+C^+C^Oj) + 


• • # 



38 


If we are using the discrete Haar transform with N = 8, 
we obtain , 

F 1 = °* 9c i + 0.1865 ( -C^+C^+C^-G^ ) 

F 2 = 0.6366(C 2 -C 3 )+0.1865(C 4 +C 5 -C 6 -C 7 ) 

F 3 = 0.6366 (C 2 +C 3 ) 

f 4 ■ 0.6366(^-05+^-07) 

If we define modified Haar coefficients as indicated in 
Fig. 3.3, these expressions get modified to, 

F x = 0.9C[ + 0.373 (C£+C£-C^-C^) 

F 2 = 0.9(C£-C£) + 0.373(C^+C£*C£-C}) (3.7) 

F 3 * 0.9(C£+CJ) 

F 4 - 0 .9 ( C^-C^+C^-C^ ) 

The rms values of the first and the second harmonic are 


given by 


*2 


•/ifTfl 


(3.8) 
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The frequency response of- the filters so obtained is 
shown in Table 3.3* 8 point Fast Haar Transformation has 

been used. The response is plotted in Fig. 3.4. The response 
is such that all undesirable frequencies will be rejected 
adequately if the operating frequency is exactly 50 Hz. With 
change in frequency, the degradation of performance is severe. 

Table 3.4 presents the simulated step response of Haar 
function filters. The results are plotted in Fig. 3.5(b). Th< 
settling time is clearly less than 20 milliseconds. Table 3.5 
and Fig. 3.5(a) show the simulated impulse response. The res- 
ponse to a centre frequency step input is shown for both fil- 
ters in Figures 3.5(c) and 3.5(d), and the simulated data in 
Tables 3.6 and 3.7 . 

3.5 Conclusion 

In this chapter, we have shown the feasibility and 
attractiveness of using Haar transform based filters for 
differential relaying. In the next one, we shall discuss the 
actual scheme that, has -been implemented. 



Table 3*3 
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Frequency Response 


Frequency 

50 Hz filter output 

100 Hz filter outp< 

(Hz). 

(per unit) 

(per unit) 

0.0 

0.0000 

0.0000 

5.0 

0.1912 

0.0948 

10.0 

0.3341 

0.1611 

15.0 

0.4010 

0.1799 

20.0 

0.4073 

0.1494 

25.0 

0.4299 

0.0932 

30.0 

0.5465 

0.0832 

35 .0 

0.7245 

0.1266 

40.0 

0.8844 

0.1415 

45.0 

0.9728 

0.0992 

46.0 

0.9799 

0.0835 

47.0 

0.9834 

0.0656 

48.0 

0.9835 

0 .0456 

49.0 

0.9803 

0.0237 

50.0 

0.9741 

0.0000 

51.0 

0.9651 

0 .0252 

52.0 

0.9537 

0.0517 

53 .0 

0.9402 

0.0793 

54.0 

0.9251 • 

0.1078 

55.0 

0.9087 

0.1367 

60.0 

0.8232 

0.2809 

65.0 

0.7626 

0.4031 

70.0 

0.7317 

0.4877 

75 .0 

0.6935 

0.5431 

80.0 

0.6112 

0.5996 

85.0 

0.4755 

0.6829 

90.0 

0.3053 

0.7842 

95.0 

0.1351 

0.8677 

96.0 

0.1042 

0.8790 

97.0 

0.0751 

0.8880 

98.0 

0.0479 

0.8946 

99.0 

0.0228 

0.8985 

LOO .00 

0.0000 

0.8999 

L01.00 

0.0204 

0.8985 

L02.0 

0.0384 

0.8946 

.03.0 

0.0537 

0.8880 

.04.0 

0 .0666 

0.8790 

.05 .0 

0.0769 

0.8677 

.10.0 

0.0969 

0.7842 

.15.0 

0.1000 

0.6829 

.20.0 

0.1436 

0.5996 

.25.0 

0.2064 

0.5431 * 

.30.0 

0.2445 

0.4877 

.35.0 

0.2358 

0.4031 

.40.0 

0.1797 

0.2809 

.45.0 

0.0923 

0.1367 

.50.0 

0.0002 

0.0000 

.55.0 

0.0705 

0.0992 
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Frequency 

50 Hz filter output 

100 Hz fi 

(Hz) 

(per unit) 

output (pe 

160.0 

0.0993 

0.1415 

165.0 

0.0816 

0.1266 

170.0 

0.0328 

0.0832 

175.0 

0.0571 

0.0932 

180 .0 

0.1167 

0.1494 

L85 .0 

0.1464 

0.1799 

L90.0 

0.1327 

0.1611 

L95.0 

0.0785 

0.0949 

>00.0 

0.0000 

0.0000 


Table 3.3 Frequency response of Haar function 
based filters 




PREQUENCY CHK) 
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Table 3.4 


Time (ms) 

Input 

50 Hz filter 

100 Hz 

0.00 

1.00 

0.17 

0.16 

2.50 

1.00 

0 .32 

0.22 

5.00 

1.00 

0.42 

0.16 

7.50 

1.00 

0.45 

0.00 

10.00 

1.00 

0.42 

0.16 

12.50 

1.00 

0.32 

0.22 

15.00 

1.00 

0.17 

0.16 

17.50 

1.00 

0.00 

0.00 

20.00 

1.00 

0.00 

0.00 

22.50 

1.00 

0.00 

0.00 

25.00 

1.00 

0.00 

0.00 

27.50 

1.00 

0.00 

0.00 

30.00 

1.00 

0.00 

0.00 


o 

© 

50.00 1.00 0.00 0.00 


Table 3.4 Step response 
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Table 3*5 


Time (ms) 

Input 

50 Hz filter 

100 Hz filte 

0.00 

1.00 

0.17 

0.16 

2.50 

0.00 

0.17 

0.16 

5.00 

0.00 

0.17 

0.16 

7.50 

0.00 

0.17 

0.16 

10.00 

0.00 

0.17 

0.16 

12.50 

0.00 

0.17 

0.16 

15.00 

0.00 

0.17 

0.16 

20.00 

0.00 

0.00 

0.00 

22.50 

0.00 

0.00 

0.00 

25.00 

0.00 

0.00 

0.00 

27.00 

0.00 

0.00 

0.00 

30.00 

• 

0.00 

0.00 

0.00 

• 

50*00 

0.00 

0.00 

0.00 


Table 3.5 Inpulse response 



Table 3<>6 


Time ( 


0.00 

2.50 
5.00 

7.50 

10.00 

12.50 

15.00 

17.50 

20.00 

22.50 

25.00 

27.50 

30.00 


Time ( 

0.00 

2.50 
5.00 

7.50 
10.00 
12.50. 

15.00 

17.50 

20.00 

22.50 

25.00 

27.50 

30.00 


1 Input 50 


0.00 

1.00 

1.41 

1. 00 

0.00 

- 1.00 

-1.41 

- 1.00 

0.00 

1.00 
1.41 
1.00 
0.00 


Hz filter 100 Hz filter 


0.00 

0.00 

0.17 

0.16 

0.39 

0.28 

0.49 

0.22 

0.49 

0.22 

0.62 

0.28 

0.86 

0.16 

0.97 

0.00 

0.97 

0.00 

0.97 

0.00 

0.97 

0.00 

0.97 

0.00 

0.97 

0.00 


Table 3.6 50 Hz sine wave step response 


Table 3.7 


2 Input 50 

0.00 

1.41 
0.00 
-1.41 
0.00 

1.41 
. 0.00 
-1.41 
0.00 

1.41 
0.00 
-1.41 
0.00 


Hz filter 100 Hz filter 


0.00 

0.00 

0.24 

0.22 

0.24 

0.22 

0.34 

0.45 

0.34 

0.45 

0.24 

0.67 

0.24 

0.67 

0.00 

0.90 

0.00 

0.90 

0.00 

0.90 

0.00 

0.90 

0.00 

0.90 

0.00 

0.90 


Table 3.7 100 Hz sine wave step response 
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CHAPTER 4 

PROPOSED SCHEME 


4.1 Introduction 

In the previous chapters* we have discussed the general 
requirements of a differential relay for transformer protection. 
Some partial solutions of the problems encountered were also 
discussed. This chapter describes the actual integrated solu- 
tion that has been implemented to overcome most of the difficul- 
ties envisaged, and presents the results of some realistic 
studies carried out to determine the performance of the proposed 
relaying scheme. It has been shown that the relay fabricated 
would function satisfactorily even in the worst possible condi- 
tions. ‘The limitations of the relay developed are also brought 
out. 

The flexibility and comprehensiveness of the proposed 
relaying scheme is limited by the abilities of the 2920 Analog 
Signal Processor. This single chip microcomputer has been 
designed especially to process real time analog signals. It 
has on-board program memory, scratchpad memory, D/A circuitry, 
A/D circuitry, digital processor, and I/O circuitry. It can 
execute programs up to 13,000 times a second if the entire 
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program memory of 192 instructions is used, which allows 
processing of signals upto 5 kHz bandwidth. With shorter 
programs, the processing speed can be increased further. 

While the presence of so many features on a single chip 
greatly helps in simplifying hardware requirements, it gene- 
rates limitations on the actual performance and ability of 
each facility provided. For example, whereas the 2920 can 
handle upto 4 analog inputs and 8 analog outputs simultaneously, 
it takes around 40 instructions for each analog input operation 
(which can be carried out side-by-side with other computations) 
and around 12 instructions for each analog output operation in 
a typical hardware configuration used. This, coupled with 
some other constraints, makes it virtually impossible to use 
all the inputs and outputs simultaneously. This prevented us 
from outputting the operating and restraining signals. 

The most severe limitation faced was of program and 
scratchpad memory limitation. The 2920 has only 192 instru- 
ctions program memory and 40 RAM locations as scratchpad memory. 
If we try to compute 16 sample Haar transforms of primary and 
secondary currents, then 32 RAM locations are required just for 
storing samples of both these currents, and it is not possible 
to perform all other computations using just 8 memory locations. 
As a result, we were limited t© 8 sample Haar transforms, which 
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are inadequate for extracting the 5th harmonic component. This 
forced us to drop over-voltage inrush detection facility in the 
scheme. Some other compromises were also necessary, and they 
will gradually become clear as we describe the implemented 
scheme step-by-step. Let us, then, move on to a description of 
what has been achieved. 

4.2 Description of the scheme 

4.2.1 Input 

The currents 1^ and I 2 » obtained from the current trans- 
formers, are first converted into voltages by passing them 
through suitable air core transactors. After this, they are 
sampled by the relay. 

The two currents are sampled within 50 microseconds of 
each other in the proposed scheme. 

4.2.2 Tap changing 

To account for the mismatch created by tap changing in a 
power transformer, one of the sampled currents is multiplied by 
a correction factor. Besides 1^ and Ig* a third analog input 

indicating the tap setting is provided to the relay. The value 

* 

of this third input TAP can be varied between +1.0 and -1.0, and 
the correction factor used is (l+TAP/8). In other wofcds, 
samples of I x are replaced by I I (l+TAP/8). With this provision. 
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tap setting changes upto +12.5?: from the nominal value can be 
accounted for. 


4.2.3 Filtering 


As mentioned in the introductory remarks, 8 sample Haar 
transforms have been used. First of all, the mean through 


I i+I 


are 


current ( — — ) , and (I^-Ig) to provide operating signals, 

obtained from I. (corrected for tap setting) and I 2 * After 

Ij+*2. 

this, these are passed through filters to obtained (— 2 * 

(ij-IgK and (12-12)2* These three signals have been designated 
restraining (RES), operating (OP) and inrush (INR) signals 
respectively. 


While implementing equations (3.7) using 2920, ii 
great advantage if the coefficients of C^, ...» C£ can be 
written as sums and differences of a small number of terms 
involving powers of 2. For instance, the binary representa- 
tion of 0.875 is 0.111, which means 0.875 = 2~ 1 +2~ 2 +2~ 5 • How ~ 
ever, if we recognize that 0.875 = 1-0.125 * l-2~ 2 , we can 
implement this coefficient in a smaller number of instructions. 
Therefore, equations (3.7) have been modified as given below t 

FI * 0.90625 q + 0.375 (C|+C£-q-q) 

F2 * 0.9023(CJ-q) + 0.375(q+C£-C£-q) ^ 

F3 » 0.9023(q+Cj) 

F4 « o.9023(q-q+c^-q) 
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Due to the change in the values of coefficients, the 
frequency response of the filters changes slightly. Table 4.1 
and Fig. 4.0 show the modified frequency response. 

4.2.4 RMS value computation 

According to equations (3.8), the rms values of the 
fundamental and second harmonic contents are given by 


A x -JA + A 

A ^ 


(3.8) 


These equations require squaring and square rooting 
operations, which require very large computational efforts. 

To avoid this, an approximation is being employed. 

/ ~~2 2 

x +y and x and y 

are both positive. Since both of these appear squared in the 
expression for 0, their polarity (sign) is not really important 
and we can use their absolute values (i.e. modulus). Conse- 
quently, we can confine out attention to the first quadrant of 
the x-y plane. 

The first quadrant is further divided into two zones by 
the line L (x * y) shown in Fig. 4.1. In region I, x 2 y * 
whereas in region II, y 2 x. 



Table 4.1 


Frequency 

50 Hz Filter 

100 Hz Filter 

(Hz) 

output (p.u.) 

output (p.u.) 


0.0 

0.0000 

0.0000 

5.0 

0.1925 

0.0951 

10.0 

0.3364 

0.1615 

15.0 

0.4037 

0.1804 

20 .0 

0.4096 

0.1498 

25.0 

0.4315 

0.0934 

30.0 

0.5482 

0.0834 

35.0 

0.7273 

0.1269 

40.0 

0.8886 

0.1419 

45.0 

0.9782 

0.0995 

46.0 

0.9854 

0.0837 

47.0 

0.9891 

0.0658 

48.0 

0.9893 

0.0457 

49.0 

0.9862 

0.0237 

50.0 

0.9801 

0.0000 

51.0 

0.9712 

0.0253 

52.0 

0.9598 

0.0519 

53.0 

0.9463 

0.0795 

54.0 

0.9311 

0.1080 

55.0 

0.9147 

0.1371 

60.0 

0.8285 

0.2816 

65.0 

0.7669 

0.4042 

70.0 

0.7348 

0.4890 

75.0 

0.6958 

0.5445 

80.0 

0.6130 

0.6012 

85 .0 

0.4769 

0.6847 

90.0 

0.3062 

0.7862 

95.0 

0.1354 

0.8700 

96.0 

0.1045 

0.8813 

97.0 

0.0753 

0.8903 

98.0 

0.0481 

0.8969 

99 .0 

0.0229 

0.9009 

-100.0 

0.0000 

0.9022 

101.0 

0.0205 

0.9009 

102.0 

0.0385 

0.8969 

•-O3.0 

0.0539 

0.8903 

1,04.0 

0.0668 

0.8813 

>-05.0 

0.0772 

0.8700 

..10.0 

0.0972 

0.7862 



115.0 

0.1004 

0.6847 

120.0 

0.1439 

0.6012 

125.0 

0.2068 

0.5445 

130 .0 

0.2447 

0.4890 

135.0 

0.3259 

0.4042 

140.0 

0.1794 

0.2816 

145.0 

0.0917 

0.1371 

150.0 

0.0011 

0.0000 

155.0 

0.0715 

0.0995 

160.0 

0.1002 

0.1419 

165.0 

0.0823 

0.1269 

170.0 

0.0333 

0.0834 

175.0 

0.0576 

0.0934 

180.0 

0.1174 

0.1498 

185.0 

0.1472 

0.1804 

190.0 

0.1335 

0.1615 

195.0 

0.0789 

0.0951 

200.0 

0.0000 

0.0000 


Table 4.1 Frequency response of Haar filters 
used (Input waveform has unity 
rms value) 




FREQUENCY CH «) -*• 

««• 4 0 < F«au€Ncr response of haar filters obtained 
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Consider an estimator of the form 0 = Ax + By in the 

region I. Without any loss of generality, we may assume that 

o 2 

(x,y) lie on the unit circle x +y = 1 (curve c) for the purpose 
of determining the accuracy of this estimator, since any other ■ 
combination of (x,y) can be obtained merely by scaling them. 

The error in the estimated value is given by, 

[ 

* (Ax + By) - 1 ( V x 2 +y 2 = 1) 

I 

At which point cm c does this error become maximum? To 

I 

I 

answer this, let us switch over to polar coordinates and 
redefine error as, 

f(0) * (A cos© + B sin©) - 1 


On differentiating with respect to 0, we obtain 
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f» (e) = -A sine + B cos© 

Therefore, f ' (0) => A sin© = B cos© => tan© = B/a 

In region I, O < 8 < 45°, or 0 <_ tan© 1. This means 
that A and B are of the same sign (positive), and B < A, if 
f(©) is to have an extremum in region I. 

The second derivative f"(©) = (-A cos© — B sin©) is non- 
positive in the region I. This means that f(©) attains a maximum 
at © = tan~^(B/A), and the maximum error is found to be, 

f(e) max = 

All these results can be summarized as given below : 

/ o o' 1- 

Objective function :©=/x+y, for x > y > 0 
Estimator used : © = AX + By, with A > B > O. 


© 

0 ° 

45° 


Error 


Comments 


A-l 

£££ - i 

r* 


.tan” 1 (B/A) 


/ A 2 +B 2 - 1 


Maximum error 
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We want to find A,B such that error is reasonably small in the 
entire region I • 


Since we are constrained to select values which are easily 
realizable , we have chosen A = 1 and B = 1/2. For these values » 
the estimator always overestimates the objective function, the 
extent of overestimation varying from OyC to 11.8^ to b.Olyi as 
0 varies from 0° to 26.56° to 45° respectively. Better estima- 
tors are, of course, possible. For example, A = 1 and B = V2-1 
(= 0.414) gives a maximum error of +8.24^ only; and A = 0.95, 

B *» 0.464 gives maximum errors of -5^ and +5.735° only. However, 
as pointed out earlier, A = 1 and B = 0.5 are easy to implement. 

Using this estimator, equations (3.8) get modified to the 
following form : 


A^ = max ( i Fjl j , | F^, t ) +5 min( i Fj^ | , | F 2 j ) 
Aj = »ax( | F 3 1 , |F 4 |) + | min( |F 3 | , |F 4 |) 


(4.2) 


These equations give rise to a maximum error of 11.8?i in 
the ratio (not 23.6 •/. because both the errors can be 

positive only) . This error has to be absorbed in the margin of 
safety in the tripping criterion. 

4.2.5 Variable bias factor 


The problems of CT mismatch etc. are less severe when the 
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currents 1^ and I ^ are small* This means that we can use a 
smaller value of bias factor S when the mean through current 
I l +I 2 

( — 5 — ) is small, leading to more sensitive relaying. However, 
S has to be increased as the mean through current becomes 
larger and the CTs begin to get saturated. 

We have used a bias factor which is directly proportional 
to the mean through current and is given by 

I i +I o 

S = (-±^) rms + 0.0625 (4.3) 

It is to be noted that for the relay developed, S can 
vary between 0.0625 and nearly 0.77 under the worst possible 
conditions. For normal load currents, S will be nearly 0.10 
to 0.72. 

4.3 Hardware description 

The hardware configuration used has been taken from the 
2920 Analog Signal Processor Design Handbook [31] and is shown 
in Fig. 4.2. 

Fig. 4.3 shows the PCB layout used, and Fig. 4.4 shows a 
photograph of the complete relay. The relay requires +5V and 
-5 V power supplies, and has no provision for input conditioning. 





Fig. 4.2 Circuit diagram of the relay 
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Fig. 4.4 Photograph of the relay which has been developed 
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4.4 Software used 

Table 4.2 gives the complete assembly language program 
listing used. The details and basic facts about 2920, as 
well as information about the assembly language, can be 
obtained from [31]. 

4.5 Performance of the proposed scheme 

The steady state tripping characteristic of a percentage 
differential relay is a circle as shown in Chapter 2 
(Section 2.3). Such a characteristic is obtained only if the 
bias factor S is a constant. For the proposed scheme, wherein 
S varies with the mean through current, any such characteristic 
cannot be obtained. We can only study a few cases and verify 
if all provisions work as stipulated. 

The other aspect that needs to be tested is the transient 
response, i.e., how soon does the relay trip after the occu- 
rrence of a fault and how do the operating, inrush, and 
restraining signals vary in the meantime. The actual hardware 
developed would provide only one output - the trip signal, 
hiding much relevant information. Therefore, we have used a 
software package SM2920, which was developed by Intel Corpora- 
tion to simulate the performance of 2920 software. The use of 
this simulator allows us to monitor signals internal to the 
relay, which would be ‘invisible* in the hardware. 
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Table 4.2 


ROM OOO = LDA DhR . CTR- ROO, NOF 
ROM 001 = SUE: DAR, KF1 . Ro2, NOF-' 
ROh 002 = LDA . CTR. KP6, ROO, CNDS 
ROM 003 = SUE DAR. DAR, ROD, I NO 
ROM 004 = LDA .CO/. NO. ROO/ INO 
ROM 005 = ADD . CO... Nl , ROO, INO 

ROM 006 = LDA . Cl,. N2, ROO, INO 

ROM 007 = ADD . Cl, . N3, ROO, INO 

ROM 008 = LDA . C6, . N4, ROO, INO 

ROM 009 = ADD . C6, . N5, ROO, NOP 

ROM 010 = LDA . C7» . N6, ROO, NOP 

ROM Oil = ADD . C7, , N7, ROO, CVTS 
ROM 012 = LDA . C4, . CO, ROO, NOP 

ROM 013 = ADD . C4, . Cl, ROO, NOP 

ROM 014 = LDA . C5, . C6, ROO, CVT7 

ROM 015 = ADD . C5, . C7, ROO, NOP 
ROM 016 = LDA . C2- . .CO, ROO, NOP 
ROM 017 = SUE C2, . Cl, ROO, CVT6 
ROM 018 = LDA . C3> . C6, ROO, NOP 
ROM C 1 9 - SUE . C3, C7, ROO, NOP 

ROM 020 = LDA . Cl, . C4, ROO, CVTS 
ROM 021 = SUE .Cl,. C5, ROO, NOP 
ROM 022 = LDA C4, . NO, ROO, NOP 
ROM 023 = SUB . C4, . Nl, ROO, CVT4 

ROM 024 = LDA . C5, . N2, ROO, NOP 

ROM 025 = SUE . C5, . N3, ROO, NOP 
ROM 026 = LDA . C6, . N4, ROO, CVTS 
ROM 027 = SUB C6, . N5. ROO, NOP 
ROM “ 028 = LDA C7, . N6, ROO, NOP 
ROM 029 = SUB . C7 , . N7, ROO, CVT2 
ROM 030 = LDA . FI, . C5, R01, NOP 
ROM 031 = ADD . FI , . C6, R01, NOP 
ROM 032 = SUB . FI, . C4, R01, CVT1 
ROM 033 = SUB . FI, . C7, R01,N0P 
ROM 034 = SUB . FI, . FI, R02, NOP 
ROM 035 = ADD .FI,. Cl, ROO, CVTO 
ROM 036 = SUB . FI, . Cl, R03, NOP 
ROM 037 = ADD . FI, . Cl, R05, NOP 
ROM 038 =' LDA . TAP, DAR, R06, NOP 
ROM 039 « SUB DAR, DAR, ROO, INI 
ROM 040 = LDA . F2, . C4, R01, INI 
ROM 041 = ADD . F2, . C5, R01, INI 
ROM 042 * SUB . F2, . C6, R01, INI 
ROM 043 = SUB . F2, . C7, R01, INI 
ROM 044 = SUB . F2, . F2, R02, INI 
ROM 045 « SUB . C2, . C2, R03, NOP 
ROM 046 * SUB . C3, . C3, R03, NOP 
ROM 047 =' ADD . C2, . C2, R05, CVTS 


coiitd.* 
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F 

046 = 

ADD 

% i 

C3-R05, NOf 

Ron 

049 - 

ADD 

^ I , 

C 1 , ROO, NOF' 

ROM 

050 = 

SUB 

F2, 

C3, ROO, CVT / 

ROM 

05 1 = 

LDA 

. F 3, 

. C2, ROO, NOP 

RDM 

052 = 

ADD 

. F3» 

. C3, ROO, NOP 

ROM 

053 = 

LDA 

. F4, 

. C4, POO, CVT 6 

ROM 

054 = 

SUB 

. F 4) 

C5, ROO, NOP 

ROM 

055 = 

ADD 

. F4, 

C6, ROO, NOP 

ROM 

056 = 

SUB 

F4» 

. C7, ROO, CVT5 

ROM 

057 = 

SUB 

. F4 f 

F4, R03, NOP 

ROM 

058 = 

ADD 

. F4.. 

. F4, R05, NOP 

ROM 

059 = 

LDA 

. CO.. 

. MO. ROO, CVT4 

ROM 

060 = 

ADD 

. CO, 

Ml, ROO, NOP 

ROM 

06 i = 

LDA 

. Q1 , 

. M2, ROO, NOP 

ROM 

062 = 

ADD 

. Cl, 

. M3, ROO, CVT3 

ROM 

063 = 

LDA 

C6, 

. M4, ROO, NOP 

ROM 

064 = 

ADD 

. C6, 

. M5, ROO, NOP 

ROM 

065 = 

LDA 

. C7, 

. M6. ROO, CVT2 

ROM 

066 = 

ADD 

. C7, 

. M7, ROO, NOP 

ROM 

067 = 

LDA 

. L4i 

. CO, ROO, NOP 

ROM 

068 = 

ADD 

. C 4 . 

. Cl. ROO, CVT1 

R JM 

069 = 

LDA 

. C5, 

. C© > ROO, NOP 

ROM 

070 = 

ADD 

. C5, 

. C7 , ROO, NUP 

ROM 

071 = 

LDA 

. C2> 

. CO, ROO, CVTO 

ROM 

072 = 

SUB 

. C2, 

.Cl, ROO, NOP 

ROM 

073 = 

LDA 

. C3 , 

. C6, ROO, NOP 

ROM 

074 = 

LDA 

. 11, 

DAR, R03, NOP 

ROM 

075 = 

ADD 

. 11, 

. TAP, R01, CND7 

ROM 

076 = 

ADD 

. 11, 

. TAP, R02, CND6 

ROM 

077 = 

ADD 

. 11, 

. TAP, R03, CND5 

ROM 

078 = 

ADD 

. 11, 

. TAP, R04, CND4 

ROM 

079 = 

ADD 

. 11, 

. TAP. R05, CND3 

ROM 

080 = 

ADD 

. 11, 

. TAP, R06, CND2 

ROM 

08 1 — 

ADD 

. 11, 

. TAP, R07, CND1 

ROM 

082 = 

ADD 

. 11, 

. TAP, R08, CNDO 

ROM 

083 = 

SUB 

. TAP 

, . TAP, L01, NOP 

ROM 

084 = 

ADD 

. 11,. 

. TAP, ROO, CNDS 

ROM 

085 * 

SUB 

DAR, ! 

DAR) ROO, IN2 

ROM 

086 = 

SUB 

. C3, . 

. C7, ROO, IN2 

ROM 

087 * 

LDA 

. Cl, . 

C4, ROO, IN2 

ROM 

088 = 

SUE! 

. Cl, . 

. C5, ROO, IN2 

ROM 

089 = 

ABS 

. OP, . 

FI, L01, IN2 

ROM 

090 » 

ABA 

. OP, . 

F2, ROO, IN2 

ROM 

091 * 

ABS 

. F2, . 

F2, ROO, NOP 

ROM 

092 » 

ABA 

. F2, . 

FI, R01, NOP 

ROM 

093 “ 

ABS 

. I NR; 

, . F3, ROO, CVTS 

ROM 

094 - 

ABA 

. INR, 

, . F4, R01, NQF 

ROM 

095 * 

ABS 

. F4, . 

F4, ROO, NOP 


contd . • 
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F OM 

U96 

=- 

AE.K 

F4> F3j R01 , CVT” 

'ROM 

097 

5= 

LDA 

. C4. MOiROliNOR 

ROM 

098 

SS 

SuL- 

C4< Ml , R01 , NOP 

ROM 

099 

sr 

SUB 

. C4» . M4, Pol, CVT 6 

ROM 

1 00 

SS 

ADD 

C4, M5. ROliNOP 

ROM 

101 

sr.: 

LDA 

. C5, . M2.. R01, NOP 

ROM 

102 

= 

SUE 

. C5 . M3. R01, CVT 5 

ROM 

103 

ss 

SUE 

. C5, . M6, R01, NOP 

ROM 

104 

ss 

ADD 

. C5, . M7.* R01, NOP 

ROM 

105 

s= 

SUE 

C4, . C4, R02, CVT 4 

ROM 

106 

ss 

SUE 

. C5, . 05, R02, NOP 

ROM 

107 

s= 

LDA 

.FI, Cl, ROD, NOP 

ROM 

108 

ss 

SUE 

. FI, . FI, R03, CVT 3 

ROM 

109 

ss 

ADD 

. FI, . FI, R05, NOP 

ROM 

110 

ss 

SUE 

. FI, . C4, ROO, NOP 

ROM 

111 

ss 

ADD 

. FI, . C5, ROO, CVT2 

ROM 

112 

ss 

LDA 

. F3, . C2, ROO, NOP 

ROM 

113 

ss 

SUE 

. F3, . C3, ROO, NOP 

ROM 

114 

ss 

SUB 

. F3, . F3, R03, CVT1 

ROM 

115 

ss 

ADD 

. F3, . F3, R05, NOP 

ROM 

116 

ss 

ADD 

. F3, . C4. ROO, NOP 

ROM 

117 

S= 

ADD 

. F-3, . L-5, ROO, UVTO 

ROM 

118 

ss 

ABS 

. AVG, . FI, ROO, NOP 

ROM 

119 

ss 

AEA 

. AVG • . F3, R01, NOP 

ROM 

120 

ss 

ABS 

. F3> . F3, ROO, NOP 

ROM 

121 

ss 

ABA 

. F3, . FI, R01, NOP 

ROM 

122 

ss 

LDA 

. TAP, . 11, R01, NOP 

ROM 

123 

ss 

ADD 

. TAP, DAR, R04, NOP 

ROM 

124 

ss 

SUB 

. 11, DAR, R03, NOP 

ROM 

125 

SS 

ADD 

. CTR, KP1, R02, NOP 

ROM 

126 

ss 

LDA 

DAR, KP6, ROO, NOP 

ROM 

127 

sr 

SUE 

DAR. DAR, R06, NUK 

ROM 

128 

ss 

SUB 

DAR, . CTR, ROO, NOP 

ROM 

129 

ss 

LDA 

. CTR, KP1, ROO, CNDS 

ROM 

130 

ss 

LDA 

. NO, . Nl, ROO, CNDS 

ROM 

131 

ss 

LDA 

. Nl, N2, ROO, CNDS 

ROM 

132 

ss 

LDA 

. N2, . N3, ROO, CNDS 

ROM 

133 

ss 

LDA 

. N3, . N4, ROO, CNDS 

ROM 

134 

ss 

LDA 

. N4, . N5, ROO, CNDS 

ROM 

1 35 

ss 

LDA 

. N5, . N6, ROO, CNDS 

ROM 

136 

ss 

LDA 

. N6, . N7, ROO, CNDS 

ROM 

137 

ss 

LDA 

. MO, . Ml, ROO, CNDS 

ROM 

138 

ss 

LDA 

. Ml, . M2, ROO, CNDS 

ROM 

139 

ss 

LDA 

. M2, . M3, ROO, CNDS 

ROM 

140 

ss 

LDA 

. M3, . M4, ROO, CNDS 

ROM 

141 

ss 

LDA 

. M4, . MS, ROO, CNDS 

ROM 

142 

2 

LDA 

. M5, . M6, ROO, CNDS 

ROM 

143 

ss 

LDA 

. M6, . M7, ROO, CNDS 


contd.; 
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ROM 

144 

=1 

LDA 

. M”. TAr.RC--. 

, CNL' : 

ROM 

145 


LDA 

. NT, . 11. ROl t 

-NDS 

ROM 

14e. 

= 

ADD 

. N7, . OLD, ROl, 

■ CNDS 

ROM 

147 

ss 

SUB 

. 11.. 11, L01 . NOP 

ROM 

148 

= 

LDA 

OLD, .11, ROO, 

CNDS 

ROM 

149 

= 

LDA 

DAR, . OP, ROO, NOP 

ROM 

150 

= 

SUB 

DAR,,. F2, L01, NOP 

ROM 

151 

= 

LDA 

. OP,'. F2, L01, CNDS 

ROM 

152 

= 

ADD 

. OP, . OP, R02, NOP 

ROM 

153 

=s 

LDA 

DAR/ . AVG> ROO/ 

NOP 

ROM 

154 


SUE 

DAR, . F3» ROO, NOP 

ROM 

i55 

=r. 

LDA 

. AVG, . F3, ROO, 

CNDS 

ROM 

156 

rs 

LDA 

DAR, . I NR, ROO, NOP 

ROM 

157 

= 

SUB 

DAR, . F4, ROO, NOP 

ROM 

158 

= 

LDA 

. I NR, . F4, ROO, 

CNDS 

ROM 

159 


ADD 

. I NR, . I NR, ROl 

, NOP 

ROM 

160 

= 

LDA 

DAR, . AVG, ROO, 

NOP 


ROM 161 = ADD DAR, KP1. ROl, NOP 
ROM 162 = ADD . RES, . AVG, ROl , CND7 

ROM 163 = ADD . RES, . AVG, R02, CNDS 

ROM 164 = ADD . RES. . AVG, R03, CND5 

ROM 165 = ADD . RES. . AVG. R04, CND4 

ROM 166 = ADD . RES. . AVG, R05, CND3 

ROM 167 = ADD . RES, . AVG, R06, CND2 

ROM 168 = ADD . RES, . AVG, R07, CND1 

ROM 169 = ADD . RES, . AVG, R08, CNDO 

ROM 170 = SUE . AVG, . AVG, L01, NOP 

ROM 171 = ADD . RES, . AVG, ROD, CNDS 

ROM 172 = SUE . OP. KP1.R05. NOP 
ROM 173 = LDA . TRIP, KP7. ROO, NOP 
ROM 174 = LDA DAR, . OP, ROl. NOP 
ROM 175 = SUE DAR,. I NR, L02, NOP 
ROM 176 = LDA . TRIP. KPO, ROO. CNDS 
ROM 177 = LDA DAR, . OP, ROO, NOP 
ROM 178 = SUB DAR, . RES, ROO, NOP 
ROM 179 = LDA . TRIP, KPO, ROO, CNDS 
ROM ISO = LDA DAR, . TRIP, ROO, NOP 

ROM 181 * LDA . TAP, . RES, ROO, NOP 

ROM 182 = SUB . RES, . RES, ROO, NOP 

ROM 183 = SUE . AVG, . AVG, L01, NOP 

ROM 184 = LDA . CTR, . CTR, ROO, NOP 

ROM 185 = LDA . CTR, . CTR, ROO, NOP 

ROM 186 - LDA . CTR, . CTR, ROO, OUTO 

ROM 187 = LDA . CTR, . CTR, ROO, OUTO 

ROM 188 = LDA . CTR, . CTR, ROO, OUTO 

ROM 189 - LDA . CTR, . CTR, ROO, OUTO 

ROM 190 * LDA . CTR,. CTR, ROO, OUTO 

ROM 191 « LDA . CTR, . CTR, ROO, OUTO 


Tabla 4*2 Software dereloped for tha relay 



67 


The transient response has been tested for three different 

i 

inputs, representing severe internal fault current, worst case 
one-phase inrush current, and worst case three-phase inrush 
current . 

Table 4.3 and Fig. 4.5 show the results for simulated fault 
conditions in an unloaded single-phase transformer. Data to 
represent samples of differential current in this situation were 
obtained from 

I x = sin (cot- hi ) s in(9- "V ) 

where -S) = arc tan Q 
Q = o)L/R 

Values of Q = 10 and 6=0 were used to obtain the fault waveform. 
It was also assumed that I 2 = 0, because the fault occurs on an 
unloaded transformer. This approach has been adapted from [11]. 

The results indicate that the tripping command was issued 
within a cycle (17.75 ms), and there were no significant variation* 
in the internal variables after this. The tripping criterion were 
satisfied with a very large margin. 

The next input considered was that of inrush current in a 
single phase transformer. Modelling of this current follows the 
approach suggested in [ 22 ] » and samples to represent i^ were 


obtained from 



Table 4*3 
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Table 4.3 Simulated results for relay operation during a fault 
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= POS (0.5 - cos A>t) 

where POS (x) = x if x 2 0, and 0 otherwise. 

This waveform corresponds to a residual flux density of 0.9 p.u. 
and a saturation flux density of 1.4 p.u. It has been assumed 
that the transformer is switched in at 0° in the voltage cycle, 
as this corresponds to the maximum inrush current. It has been 
verified [22] that this waveform has the minimum ontainable amount 
of second harmonic (17.1j<) for any residual flux level. I 2 was 
assumed to be zero. 

The results are given in Table 4.4 and Fig. 4.6. Again, we 
find that all internal variables stabilize within a cycle, and ! 

i 

tripping is avoided by a comfortable margin for this worst case 
single phase inrush current. 

The worst inrush current waveform for a three phase trans- 

| 

former (for the purpose of relaying) has been taken to be 

i 

i x = POS (0.5 - cos wt) + POS (sin («t-30°)) 
and i~ = 0 . J 

I 

This has, again, been taken from [22], and has, a second harmonic f 
content of only 16.5?i. The results, given in Table 4.5 and j 
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Fig. 4.7, again show that tripping is avoided by a comfortable 
margin. 

Table 4.6 shows the results of some studies designed to 
test various new features of the proposed scheme. A brief look 
at the last two columns shows that the variable bias feature is 
working, exactly as envisaged. Another point to be noted is 
that the relay always overestimates the rms value of the mean 
through current, as asserted in Section 4.2.4. 

In the first simulation, tripping occurs as expected. The 
next simulation is similar in nature, and was designed to test 
the variable bias factor provision. In the third simulation, 
the TAP input is set to -0.8. This reduces Ij_ by 10/. f and 
tripping does not occur. Under the same conditions, a phase 
difference is introduced, in the fourth simulation, between 1^ 
and I 2 , which again causes tripping. 

The first three results of table 4.6 were verified on the 
hardware that has been developed. 

4.6 Conclusion 

We have found that all features described in the earlier 
part of this chapter work satisfactorily. The relay performance 
is expected to deteriorate sharply as the operating frequency 



Table 4.D 
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varies from 50 Hz, although this has not been verified 
experimentally. 

The results obtained indicate that the relay can be 
expected to be highly reliable and fast in operation. 
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CHAPTER 5 

CONCLUSION 


A percentage differential relay for transformer protection 
using harmonic restraint to avoid tripping due to magnetizing 
inrush current has been designed, fabricated, and tested. The 
hardware realization uses the Intel 2920 Analog Signal Proce- 
ssor. Digital signal processing techniques have been used, 
giving more flexibility, reliability, and higher speed of 
operation. 

The main objective of hardware simplification has been 
achieved to a very large extent. This will lead to ease in 
maintenance and repair. Identical units have been used for all 
the phases of a three-phase transformer, which would again help 
in maintenance. The protection of each phase is independent of 
the other phases. 

Several unique features have been added to the conventional 
scheme of unit protection. These include provision for taking 
tap changing into account, and variable bias factor for greater 
sensitivity. A new method of approximating the rms value of 
harmonics from the corresponding sine and cosine wave amplitudes 
has been tried out successfully. 
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The relay requires very little power to operate, and 
would impose very low burden on the current transformers used 
for monitoring winding currents small changes in component 
values are hardly expected to cause any change in the perfor- 
mance of the relay (i.e., the relay is insensitive to component 
values). Simple hardware would lead to low cost also. 

The tests and simulations performed reveal that the relay 
can be expected to perform satisfactorily in actual operation. 

One of the major limitations of the proposed relay is its 
high sensitivity to frequency fluctuations. A careful look at 
the frequency, response of Haar filters (Fig. 3.4) reveals that 
the attenuation of undesirable harmonics decreases sharply as 
the operating frequency changes from 50 Hz. As a result, the 
relay may not perform as envisaged if the operating frequency 
changes significantly. 

Some provision to make the relay self checking could also 
be introduced. 

Another drawback is the use of 8 sample Haar transforms 
for filtering. If the phenomenon of over voltage inrush 
(caused by switching) is also to be taken into account, we must 
use 16 sample Haar transforms. This might be possible using 
some more recent integrated circuits like Intel 8748 (single 
chip microcomputer) and some others. This possibility needs to 
be explored. 
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